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ABSTRACT. We study heat semigroups generated by self-adjoint Laplace operators on 
metric graphs characterized by the property that the local scattering matrices associated 
with each vertex of the graph are independent from the spectral parameter. For such op- 
erators we prove a representation for the heat kernel as a sum over all walks with given 
initial and terminal edges. Using this representation a trace formula for heat semigroups is 
proven. Applications of the trace formula to inverse spectral and scattering problems are 
also discussed. 



1. Introduction 

Metric graphs or networks are one-dimensional piecewise linear spaces with singular- 
ities at the vertices. Alternatively, a metric graph is a metric space which can be written 
as a union of finitely many intervals, which are either compact or [0, oo); any two of these 
intervals are either disjoint or intersect only in one or both of their endpoints. It is natural 
to call the metric graph compact if all its edges have finite length. 

The increasing interest in the theory of differential operators on metric graphs is moti- 
vated mainly by two reasons. The first reason is that such operators arise in a variety of 
applications. We refer the reader to the review [34|, where a number of models arising 
in physics, chemistry, and engineering is discussed. The second reason is purely mathe- 
matical: It is intriguing to study the interrelation between the spectra of these operators 
and topological or combinatorial properties of the underlying graph. Similar interrelations 
are studied in spectral geometry for differential operators on Riemannian manifolds (see, 
e.g. [7 1, [ 15 1) and in spectral graph theory for difference operators on combinatorial graphs 
(see, e.g. J9))- Metric graphs take an intermediate position between manifolds and combi- 
natorial graphs. 

In the present work we continue the study of heat semigroups on metric graphs initiated 
in ||281 . There we provided sufficient conditions for a self-adjoint Laplace operator to 
generate a contractive semigroup. Moreover, we proved a criterion guaranteeing that this 
semigroup is positivity preserving. For earlier work on heat semigroups generated by 
Laplace operators on metric graphs and their application to spectral analysis we refer to 

na, ma, m, ma, gn, ma, ma. 

In this article we study heat semigroups generated by self-adjoint Laplace operators 
which are characterized by the property that the local scattering matrices associated with 
each vertex of the graph are independent of the spectral parameter. All boundary conditions 
leading to such operators are described in Proposition l2.4l below. In particular, Neumann, 
Dirichlet, and the so called standard boundary conditions are in this class. 
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Our main technical tool to study heat semigroups on metric graphs are walks on edges 
of the graph, a concept developed in ll28l . ll29ll . We will revisit this concept in Section [5] 
below. Furthermore, we will provide a representation for the heat kernel as a sum over all 
walks with given initial and terminal edges. This representation relates the topology of the 
graph to analytic properties of the heat semigroup. 

In Section[4]we prove a trace formula for heat semigroups on arbitrary (compact as well 
as noncompact) metric graphs, an analog of the celebrated Selberg formula for differential 
operators on Riemannian manifolds (see [38 1, |48| for the case of compact manifolds and 
ll22l . Il37l for the noncompact case). A discrete analog of the Selberg trace formula on 
fc-regular trees is discussed in [49 1. 

The trace formula expresses the trace of the semigroup difference as the sum over all 
cycles on the graph, that is, equivalence classes of closed walks. In the particular case 
of compact graphs and standard boundary conditions our result recovers the well-known 
trace formula obtained by Roth fl46), 07). Related results can be found in J35], EO], BP . 
Il50l . In the physical literature trace formulas for Laplace operators on metric graphs have 
been discussed in 0, ll3D . Il32l . Il33l . Their applications to quantum chaos and spectral 
statistics are reviewed in the recent article lfl6l . 

As an application of the trace formula, in Section|5]we discuss inverse spectral and scat- 
tering problems. The inverse problems considered here consist of determining the graph 
and its metric structure (i.e. the lengths of its edges) from the spectrum of the Laplace op- 
erator and the scattering phase (that is, half the phase of the determinant of the scattering 
matrix), under the condition that the boundary conditions at all vertices of the graph are 
supposed to be known. Another kind of the inverse scattering problem, the reconstruc- 
tion of the graph and the boundary conditions from the scattering matrix, has been solved 
recently in ||29l . 

The results of Section |5]provide a mathematically rigorous solution of the inverse scat- 
tering problem as proposed by Gutkin and Smilansky in lfl9l . Also these result extend 
the solution of the inverse spectral problem on compact graphs given by Kurasov and 
Nowaszyk in ll35l to more general boundary conditions. 

Acknowledgments. It is a pleasure to thank the organizers of the conference "Transport 
and Spectral Problems in Quantum Mechanics" held at the University of Cergy-Pontoise 
in September 2006 for a very interesting and enjoyable meeting, both scientifically and 
socially. The authors would like to thank M. Karowski for helpful discussions. 

2. Background 

A finite graph is a 4-tuple Q = (V,X, £, d), where V is a finite set of vertices, X is a 
finite set of internal edges, £ is a finite set of external edges. Elements in X U £ are called 
edges. The map d assigns to each internal edge i e Ian ordered pair of (possibly equal) 
vertices d(i) := {v\, V2] and to each external edge e G £ a single vertex v. The vertices 
vi =: d~(i) and v-i —: d + (i) are called the initial and terminal vertex of the internal 
edge i, respectively. The vertex v — d(e) is the initial vertex of the external edge e. If 
d(i) = {v, v}, that is, d~(i) = d + (i) then i is called a tadpole. A graph is called compact 
if £ = 0, otherwise it is noncompact. 

Two vertices v and v' are called adjacent if there is an internal edge i E X such that 
v G d(i) and v' G d(i). A vertex v and the (internal or external) edge j G X U £ are 
incident if v G d(j). 

We do not require the map d to be injective. In particular, any two vertices are allowed 
to be adjacent to more than one internal edge and two different external edges may be 
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incident with the same vertex. If d is injective and d~{i) ^ d + (i) for all i El, the graph 
Q is called simple. 

The degree deg(u) of the vertex v is defined as 

deg(u) = |{e G £ | 9(e) = v}\ + \{i E 1 | 0-(t) = v}\ + \{i El \ d+(i) = v}\, 

that is, it is the number of (internal or external) edges incident with the given vertex v by 
which every tadpole is counted twice. 

It is easy to extend the First Theorem of Graph Theory (see, e.g. [ 1 1 1) to the case of 
noncompact graphs: 

(2.1) Y,deg(v) = \S\+2\l\. 

vev 

A vertex is called a boundary vertex if it is incident with some external edge. The set of 
all boundary vertices will be denoted by dV. The vertices not in dV, that is in V \ dV are 
called internal vertices. 

The compact graph Q- lnt = (V,l, 0, d\j) will be called the interior of the graph Q = 
(V, 1, £ , d). It is obtained from Q by eliminating the external edges. 

The star S[y) C £ U I of the vertex v E V is the set of the edges adjacent to v. 

Throughout the whole work we will assume that the graph Q is connected, that is, for 
any v, v' € V there is an ordered sequence of vertices {vi = v, i>2, ■ ■ ■ , v n — \,v n = v'} 
such that any two successive vertices in this sequence are adjacent. In particular, this 
assumption implies that any vertex of the graph Q has nonzero degree, i.e., for any vertex 
there is at least one edge with which it is incident. 

We will endow the graph with the following metric structure. Any internal edge i e I 
will be associated with an interval [0, Oj] with <jj > such that the initial vertex of i 
corresponds to x — and the terminal one to x = a. L . Any external edge e S £ will 
be associated with a semifine [0, +oo). We call the number a, the length of the internal 
edge i. The set of lengths {aj}, e x, which will also be treated as an element of IRl x l, will 
be denoted by a. A compact or noncompact graph Q endowed with a metric structure is 
called a metric graph and is written as (G,a). 

Given a finite graph Q = (V,T, £, d) with a metric structure a = {ai}iei consider the 
Hilbert space 

(2.2) H = H(£,I,a)=H s ®H I , H £ =0H e , H T = 0Wi, 
where Hj = L 2 (Ij) with 



7 = f[o,aj-] if jex, 
3 l[0,oo) if je£. 



Let Ij be the interior of Ij, that is, Ij = (0, aj) if j el and Ij = (0, oo) if j G £. 

In the sequel the letters x and y will denote arbitrary elements of the product set 

x L. 

jeeui 

By T>j with j G £ U X denote the set of all ipj G Hj such that ipj (x) and its derivative 
ip'j(x) are absolutely continuous and ip'-{x) is square integrable. Let T>j denote the set of 
those elements ipj G T>j which satisfy 

^" (0) = for ie£ and ^(0) = ^K) = 

^(o) = o tor jet and ^(o) = ^K) = o tor jel - 
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Let A be the differential operator 

(2.3) (AV)j (as) = ^j(z), jelUf 
with domain 

V° = X>o c W. 

jefux 

It is straightforward to verify that A is a closed symmetric operator with deficiency indices 
equal to \S \ +2\1\. 

We introduce an auxiliary finite-dimensional Hilbert space 

(2.4) K = K(8,X) = fC £ ® /df } ffi /C^ +) 

with /C f Cl f I and /C^ ±} S Cl z l. Let d /C denote the "double" of /C, that is, d JC = /C © /C. 
For any E V := Pj we set 

(2.5) :=^©^' G d /C, 
with r/> and i// defined by 

K(0)} e€ £ 
W(0)}i € i 





/{Ve(0)}eeA 




(2.6) 


i= {^(o)}t€i , 






\{V ) i( a i)}iez/ 





Let J be the canonical symplectic matrix on d fC 

I 



(2.7) J 







with I being the identity operator on JC. Consider the non-degenerate Hermitian symplectic 
form 

(2.8) ^([0],[V>]):=([0],JM>, 

where (•, •) denotes the scalar product in d K = <C 2 (\ £ \+ 2 \ X \\ 

A linear subspace A4 of d JC is called isotropic if the form co vanishes identically on M.. 
An isotropic subspace is called maximal if it is not a proper subspace of a larger isotropic 
subspace. Every maximal isotropic subspace has complex dimension equal to \£ \ + 2\T\. 

Let A and B be linear maps of K. onto itself. By (A, B) we denote the linear map from 
d K, = JC K. to K. defined by the relation 

(A,B) (xi®X2) :=A X i+B X 2, 

where xi, X2 & IC. Set 

(2.9) M(A,B) := Ker(A,B). 

Theorem 2.1 ((24)). A subspace M. C d JC is maximal isotropic if and only if there exist 
linear maps A, B : fC — » K such that M. = M(A, B) and 

(i) the map (A, B) : d IC — > JC has maximal rank equal to \£ I + 2\X\, 

(2.10) . . 

(ii) AB is self-adjoint, AB^ = BA*. 

Under the conditions ( |2.10t both A ± ikB are invertible for all k > 0. 
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Definition 2.2. Two boundary conditions (A, B) and (A' , B') satisfying ( 12.1 Ob are equiv- 
alent if the corresponding maximal isotropic subspaces coincide, that is, M(A, B) = 
M(A',B'). 

The boundary conditions (A,B) and (A',B') satisfying ( 12.101 are equivalent if and 
only if there is an invertible map C : K. — > JC such that A' = CA and B' = CB (see 
Proposition 3.6 in |29l ). 

By Lemma 3.3 in 11291 . a subspace Ai (A, B) C d JC is maximal isotropic if and only if 

(2.11) M(A,B) ± = M(B,—A). 
We mention also the equalities 

M(A,B) ± = [Ker (A, B)] ± = Raa(A,B)\ 

M(A, B) = Ran(-S, J 4) t . 

There is an alternative parametrization of maximal isotropic subspaces of d JC by unitary 
transformations in JC (see [28| and Proposition 3.6 in |29|). A subspace A4(A, B) C d JC 
is maximal isotropic if and only if for an arbitrary k e M \ {0} the operator A + ikB is 
invertible and 

(2.12) 6(k; A, B) := -{A + ikS)" 1 ^ - ikB) 

is unitary. Moreover, given any k e M\{0}the correspondence between maximal isotropic 
subspaces M. C d JC and unitary operators 6(k; A, B) g U(|f | + on JC is one-to-one, 
a result dating back to Bott [6| and rediscovered in (4), ETll . and |25|. Therefore, we will 
use the notation 6(k; M) for 6(k; A, B) with M(A, B) = M. 

Under the duality transformation A4 i— > M. , as a direct consequence of ( 12.111 ) and 
( 12.12b . the operators (12.12b transform as follows (see Corollary 2.2 in 1241 ): 

(2.13) 6(k;7W ± ) = -6(k- 1 ;X). 

There is a one-to-one correspondence between all self-adjoint extensions of A° and 
maximal isotropic subspaces of d JC (see ll24l . 11291 ). In explicit terms, any self-adjoint 
extension of A is the differential operator defined by (12.3b with domain 

(2.14) Dom(A) = {ip € T>\ [ip] € M}, 

where A4 is a maximal isotropic subspace of d IC. Conversely, any maximal isotropic 
subspace M. of d K. defines through (12.141 ) a self-adjoint operator A(M,a). If X = 0, 
we will simply write A(A4). In the sequel we will call the operator A(A4,a) a Laplace 
operator on the metric graph (G,a). From the discussion above it follows immediately 
that any self-adjoint Laplace operator on Ti equals A(A4,a) for some maximal isotropic 
subspace M. Moreover, A(M,a) = A(M',a) if and only if M = M'. 

From Theorem l2.1l it follows that the domain of the Laplace operator A(A4 , a) consists 
of functions ip £ T> satisfying the boundary conditions 

(2.15) AiP_ + BtJ/ = 0, 

with (A, B) subject to d2T9l and ( 12.10b . Here ip and ip' are defined by d2~6l . 

With respect to the orthogonal decomposition fC = Kg fC^ ' K.^ any element \ 
of JC can be represented as a vector 

(2.16) X = 
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Consider the orthogonal decomposition 

(2.17) AC = 0£„ 

veV 

with C v the linear subspace of dimension deg(w) spanned by those elements ( 12.16b of JC 
which satisfy 

Xe = if e e £ is not incident with the vertex v, 

(2.18) x\ ^ = if u is not an initial vertex of i € I, 

=0 if i? is not a terminal vertex of i el. 

Obviously, the subspaces £ Ul and C V2 are orthogonal if v% ^ v-2- 

Set d £„ := C v £„ = C 2dc6 ^ v \ Obviously, each d C v inherits a symplectic structure 
from d JC in a canonical way, such that the orthogonal decomposition 

holds. 

Definition 2.3. Given the graph Q = Q(V, X, £ , 0), boundary conditions (A, B) satisfying 
(12.10b are called local on (J /ff/ze maximal isotropic subspace Ai(A, B) of d JC has an 
orthogonal symplectic decomposition 

(2.19) M(i,5) = 0M„, 

with M v maximal isotropic subspaces of d C v . Otherwise the boundary conditions are 
called non-local. 

By Proposition 4.2 in 11291 , given the graph Q — Q(V,T,£,d), the boundary conditions 
(A, B) satisfying ( 12.101 ) are local on Q if and only if there is an invertible map C : JC — > JC 
and linear transformations A(v) and B(v) in C v such that the simultaneous orthogonal 
decompositions 

(2.20) CA = A(v) and CB = 

ugy vev 
are valid. From the equality M{A, B) = M{CA, CB) it follows that the subspaces M v 
in ( 12.191 1 are equal to M(A(v), B(v)). 

Boundary conditions (A(v), B(v)) induce local boundary conditions (A, B) on the 
graph Q with 

(2.21) -4 = A ( v ) and B = 
From (12.20l i we get that 

(2.22) 6(k; A, B) = 6(k; CA, CB) = 6(k; B(v)) 

vev 

holds with respect to the orthogonal decomposition (12.17l >. 
The following proposition is taken from [29|. 

Proposition 2.4. Let M. = A4(A,B) be a maximal isotropic subspace. The following 
conditions are equivalent: 

(i) S(k; A4) is k-independent, 
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(ii) 6(k;.M) is self-adjoint for some k > 0, 

(iii) for some k > there is an orthogonal projection P such that 6(k; M) = I — IP, 

(iv) ABi = 0. 

Since this proposition will be crucial in what follows, we recall the 

Proof, (i) <^> (ii). Assume that 6(k; M) is k-independent. Then, by (12.121 1. for any eigen- 
vector x £ JC with eigenvalue A the equality 

(A + l)ylx + ik(A-l)5x = 

holds for all k > 0. Under the conditions (12.101 > we have Ker A _L Ker B (see Lemma 3.4 
in (29)). Hence, A G {—1, 1}. Thus, 6(k; M) is self-adjoint for all k > 0. 

Conversely, assume that 6(k; A4) is self-adjoint for some ko > 0. Due to the obvious 
equality 

(2.23) 6(k;M) = ((k-ko)6(ko;M) + (k + k )) _1 ((k + ko)6(ko;X) + (k-ko)), 

it is self-adjoint for all k > 0. Let \ £ JC be an arbitrary eigenvector of &(ko;A4) 
corresponding to the eigenvalue A € { — 1,1}. Observing that 

(k + k )A + k-k = 
(k-k )A + k + k 

again by ( |2.23b , we conclude that \ is an eigenvector of 6(k; M) corresponding to the 
same eigenvalue A for all k > 0. Thus, 6(k; M) does not depend on k > 0. 

The equivalence (ii) <^> (iii) is obvious. 

The equivalence (iv) <£4> (ii) follows directly from the identity 

6(k;.M) -6(k;A^) f 

= 2ik(A + ikB)- 1 [B(A* - ikB^) + (A + ikB)B i ] (A^ - ikB^y 1 
= 4ik(A + ikB^AB^A^ - ikS^ -1 . 

□ 

We will write 6(M) instead of 6(k; Ai), whenever any of the equivalent conditions of 
Proposition ^. 4| is met. Analogously we will drop the k-dependence in ( 12.221 ): 

= &(A(v),B(v)) - &{M V ). 
vev vev 

From Proposition 3.5 in ll28l it follows that for any maximal isotropic subspace M. 
satisfying any of the conditions of Proposition 12.41 the Laplace operator — A(vW,a) is 
nonnegative. 

Remark 2.5. Assume that the maximal isotropic subspace M. C d JC satisfies any of 
the conditions of Proposition \2.4\ By (12.1 11 1 the orthogonal maximal isotropic subspace 
A4 1 - C d JC then also satisfies the conditions of Proposition^^ From ( 12.131 l it follows that 
6(M ± ) = -e(M). 

Obviously, Dirichlet A = I, B = and Neumann A = 0, B = I boundary conditions 
satisfy the conditions of Proposition ^. 41 with ©(I, 0) = —I and ©(0, 1) = I, respectively. 
We now provide two important examples of boundary conditions satisfying the conditions 
referred to in Proposition l2.4l 
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Example 2.6 (Standard boundary conditions). Given a graph Q — Q(V,T, £, d) for each 
vertex v € V with deg(u) > 2 define the boundary conditions (A(v), B(v)) the deg(i>) x 
deg(u) matrices 



A(v) 



Clearly, A(v)B(v)^ — and (A(v), B(v)) has maximal rank. The corresponding unitary 



(I 


-1 


.. 





o\ 




f° 





.. 





°\ 





1 


-1 . . 
















.. 














1 . . 
















.. 


















, B(v) = 


















.. 


1 


-1 










.. 













.. 









V 


1 


1 . . 


1 


1/ 



matrices (12.12b are given by 



(2.24) 



[&(A(v),B(v))] e , e , 



deg{v) 

with 5 e . e > Kronecker symbol. If deg(w) = 1, we set A(v) — 0, B(v) = 1 (Neumann 
boundary conditions) such that ( |2.24| l remains valid. 

The local boundary conditions (A, B) on the graph Q defined by ( 12. 2U will be called 
standard boundary conditions. We use the notation M. s t for the corresponding maximal 
isotropic subspace. 

Remark 2.7. Consider a graph with no internal lines Q — ({v}, 0,£,d) and \£\ > 2. By 
Proposition 2.1 in 11211 . the set of all isotropic subspaces satisfying any of the equivalent 
conditions of Proposition \2.4\ contains precisely four spaces, which correspond to the 
boundary conditions invariant with respect to permutations of edges: M(J, 0) (Dirichlet), 
Ai(0,T) (Neumann), standard Ai s t, and co-standard M^ t . 

Furthermore, by a result in 1301 . in the set of all isotropic subspaces satisfying any of 
the equivalent conditions of Proposition |2~?1 M s t is the only one with the property that 
every function in the domain of A(Ai) is continuous at the vertex v. 

Example 2.8 (Magnetic perturbations of standard boundary conditions). If the maximal 
isotropic subspace A4(A, B) satisfies any of the equivalent conditions of Proposition f2~4\ 
then for any unitary U we have 

AU(BU)^ = AB ] = 0. 

Thus, the maximal isotropic subspace A4 U := A4(AU, BU) also satisfies the conditions 
of Proposition \2.4\ In particular, since 

&{M U ) = U ] &{M)U, 

we have the relation 

(2.25) tr K &{M u ) =tT K e(M). 

A special choice of unitary matrices U corresponds to magnetic perturbations of the 
Laplace operator A(A4, a). By a result in 11271 any magnetic perturbation of the Laplace 
operator A(Ai,a) is unitarily equivalent to A(A4 U ,a) with some U = U v , where 

every U v is unitary and diagonal with respect to the canonical basis in C v , 

U v = diag ({e^W} jeSW ) . 

In particular, any magnetic perturbation of standard boundary conditions ( see Example 
satisfies the conditions of Proposition \2~4\ 
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3. Heat Kernel and Walks on the Graph 

3.1. The Resolvent. The structure of the underlying Hilbert space H (12.21 gives naturally 
rise to the following definition of integral operators. 

Definition 3.1. The operator K on the Hilbert space Tt is called integral operator if for all 

j,j' G £ Ul there are measurable functions Kj ■) : Ij x Ij, — > C with the following 
properties 

(i) Kjjtfxj, -)ipjr(-) G L 1 (I ji) for almost all Xj G Ij, 

(ii) ijj = Kip with 

(3.1) *Pj(zj)= / K ],j'i x hVj') { Pr(yj') d V3'- 

The (\X\ + \£\) x (|X| + \£\) matrix-valued function (x,y) <— > K(x,y) with 

is called the integral kernel of the operator K. 

Below we will use the following shorthand notation for ( 13.1b : 

ip(x) = J K(x,y)ip(y)dy. 

We denote 



(3.2) i?(k;a) := I I 



, e~ ik " 

and 

/0 

(3.3) T(k;a) := e ik ^ 

\0 e ik ^ 

with respect to the orthogonal decomposition ( 12.4l i. The diagonal \X\ x \X\ matrices e ±lk - 
are given by 

(3.4) [e ±Sk& ] j k = S jk e ±Ska i for j, k G 1. 
Lemma 3.2. For any maximal isotropic subspace Ai C d K. the resolvent 

{-A(M;a) -k 2 )- 1 for k 2 G C \ spec(-A(7W; a)) with det(A + ikB) ^ 

is the integral operator with the + \£\) x (\X\ + \£\) matrix-valued integral kernel 
rM(x, y, k, a), Im k > 0, admitting the representation 

r M (x,y;k,a) = A 0) (x,y,k) 

(3.5) ; 

+ ^$0, ^Ri^a)- 1 ^ - 6(k;M)T(k;a)}- 1 &{k;M)R(k;a)- 1 <i>{y, k) T , 
where R(k] a) is defined in ( 13.21 ), the matrix Q(x, k) is given by 

<S>(xk)-=(^ X ^ ° ° 
with diagonal matrices <fi(x, k) = diagje 1 ^ }j^e > <P±{%> k) = diag{e ±lk:Ej '}jex, and 

QM x i-Vj\ 

[r [ ] {x,y, k)]jj, = i5j tj > — , x 3 , Vj G Ij . 
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If I = 0, this representation simplifies to 

r M (x, y, k) = r<°> (a:, y, k) + ±<f>(x, k)6(k; M)4>(y, k). 

The integral kernel r» (x, j/; k, a) is called Green 's function or Green 's matrix. 
The proof of Lemma |3~2l is given in |28l . 

3.2. Walks on Graphs and Cycles. We recall the following definitions from |28l . A 
nontrivial walk w on the graph Q from j G £ U X to j' 6 £ U X is an ordered sequence 
formed out of edges and vertices 

(3-6) {j> , jl,fl, • ■ -,jn,Vn,j'} 

such that 

(i) jl, ■ ■■ ,jn e T, 

(ii) the vertices vo £ V and v n £ V satisfy vq G vq G u„ € <9(j'), and 
"n G 3(j„); 

(iii) for any fc G {1, . . . , n — 1} the vertex Wfc G V satisfies Wfc £ d(jk) and Ufc S 

<9(jfc+i); 

(iv) Ufe = Ufc+i for some k £ {0, . . . , n — 1} if and only if is a tadpole. 

If j, j' £ £ this definition is equivalent to that given in [29]. 

The number n is the combinatorial length |w| com b and the number 

n 
fc=l 

is the metric length of the walk w. 

A trivial walk on the graph Q from j £ £ U X to j' £ £ U X is a triple {j, u, j'} such 
that i? G d(j) and w G Otherwise the walk is called nontrivial. In particular, if 

d{j) = {w ,wi}, then {j,v , j} and {j,vi,j} are trivial walks, whereas {j,v ,j,Vi,j} 
and {j, v% , j, Vq, j} are nontrivial walks of combinatorial length 1 and of metric length dj. 
Both the combinatorial and metric length of a trivial walk are zero. 

We will say that the walk ( 13.61 ) leaves the edge j through the vertex vq and enters the 
edge j' through the vertex v n . A trivial walk {j, v,j'} leaves j and enters j' through the 
same vertex v. 

For any given walk w from j £ £ U X to j' £ £Ulwe denote by u_ (w) the vertex 
through which the walk leaves the edge j and by v+ (w) the vertex through which the walk 
enters the edge j'. For trivial walks one has u_(w) = v+(w). 

Assume that the edges j, j' £ £ U X are not tadpoles. For a walk w from j to j' we set 

dist(:Tj, u_(w)) 

and 

dist(xj' , f+(w)) 

A walk w = {j.vo, ji,vi, . . . ,j n ,v n ,j'} traverses an internal edge i £ X if jk — i for 
some 1 < As < ri. It visits the vertex v if Vk = v for some < k < n. The score n(w) 




Xj if u_(w) =d (j), 

a,j - xj if u_(w) = 

if v+(w) =d~(j), 

dj-Xj if u+(w)=9 + (j). 
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of a walk w is the set {rii(w)}i e x with rti(w) > being the number of times the walk w 
traverses the internal edge i £ X. In particular, 

|w| = y^o f n,(w). 

We say that the walk is transmitted at the vertex Vk if either Vk = d(e) or Vk — d(e') 
or Vk £ d(ik), Vk £ d(ik+i), and ik ^ ik+i- We say that a trivial walk from e' to e is 
transmitted at the vertex v = 9(e) = d(e') if e ^ e'. Otherwise the walk is said to be 
reflected. 

Let Wj.ji, j, j' £ £ U X be the (infinite if X ^ 0) set of all walks w on Q from j to f. 
By fe)> ^ ^ (No)' 1 ' we denote set of all walks won? from j to j' with score n. 
A walk 

(3.7) W = {j,V 0> jl,Vl, . . .,j n ,v n ,j'} 

is called closed if j = j' . It is called properly closed if it is closed and v_ (w) = v ^ v n = 
v + (w). For any closed walk w we denote by j(w) its initial edge, that is, j(w) = j = j'. 

For instance, let j be an arbitrary internal edge with d(j) — {vo, Vi}, vo ^ v\. Then, 
the walk {j,Vo,j} is not properly closed, whereas {j,VQ, j,vi, j} is. Any closed walk 
from an external edge is not properly closed. 

We will say that two properly closed walks w and w' are equivalent, if they can be 
obtained from each other by successive application of the transformation of the form 

{j,V0,jl,1>l, •••,in,«n,i} -> {jl,Vl,.. ■ , jn,Vn,j,V ,jl}. 

A cycle is an equivalence class of properly closed walks. We will say that the cycle c is 
associated with a walk w and write c(w), if w is in the equivalence class c. 
The number 

(3.8) \c\:=\w\+ aj , 

where w is an arbitrary walk in the equivalence class c and j = j(w), will be called the 
metric length of the cycle c. Obviously, this definition does not depend on the particular 
choice of the walk w in c. The set of all cycles on the graph Q will be denoted by £. 

We call a cycle c primitive if for any w in c there is no integer p > 2 such that 
{p _1 ni(w)}i S i is a score of a properly closed walk. For instance, if j £ X, d(j) = 
{vo,Vi}, vq 7^ i>i, the cycle associated with the properly closed walk {j,Vo,j,vi,j} is 
primitive, whereas the cycle associated with the properly closed walk {j, Vo,j, V\, j, vq , j, 
Vi,j} is not. 

For an arbitrary cycle c and any p £ N we denote by pc the unique cycle with the follow- 
ing property: For any walk w in pc there is a walk w' in c with the score {p _1 ni(w)}j £ x- 
The set of all primitive cycles on the graph Q will be denoted by € pT - lm . 

The reverse of the walk w is the walk w rev is {j',v n ,j n , . . . , Vo,j}- It may 

happen that w ICV = w. If w is a properly closed walk, then its reverse w rcv is also 
properly closed. We will write c lov for the equivalence class associated with w rev for w 
in c. Obviously, the map c h- > c rcv satisfies (pc) xev = pc rcv for any p £ N. From what has 
been just said, it follows that the case c = c rov may occur. 

3.3. Combinatorial Expansion of the Resolvent. For any k e C with Imk > 0, the 
operator T(k; a) defined in (13.31 l is a uniform contraction. Therefore, 

oo 

(3.9) [I-6(k;A^)T(k;a)]- 1 S(k;A^) = ^ 6(k; M) (T(k; o)6(k; M)) n 

n=0 
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converges uniformly in k for k in any of the sets {k £ C|Im k > e > 0}. Inserting ( 13. 9\ 
into ( 13.5b . we get 

r M {x,y;k,a) = r (0) (:r,?/, k) 

OO 

+ ^ S ^ ^) (r(k;a)6(k; M)) n ^(k;^)- 1 *^, k) T . 

For maximal isotropic subspaces satisfying any of the equivalent conditions of Proposition 
31 ©(k; A4) is independent of k. Thus, using (13.91 ) we get 



Proposition 3.3. Assume that the graph Q has no tadpoles. For any maximal isotropic sub- 
space M. satisfying any of the equivalent conditions of Proposition^^ the Green function 
of the Laplace operator A(M,a) has the absolutely converging expansion 



[r M (x,y;k,a)] jtjl = ^ s /^ ' Vjl 

(3 ' 10) + A_ e ikdist( ^^- (w)) VKM(w)e ik l w le ikdist( ^'' t,+ (w)) , lmk>0, 

2k 



where Wm ( w ) is a (complex-valued) weight of the walkw = {j,VQ,ji,Vi, . . . ,jn,v n ,j'}, 

|w| comb 

(3.ii) w M (w)= [J [e(^(«0>-B(«,))] <l+lA . 

3 .4. Heat Kernel. The semigroup generated by the positive operator — A (M , a) is related 
to its resolvent by the Dunford-Taylor integral (see ||23l Section IX. 1.6]) 

1 



27ri 



c-- ■ - =-— I e tx i A(M,a) A) d\. 



7 



where 7 is any contour encircling a positive semiline counterclockwise. The integral con- 
verges in the sense of Bochner. Using the well-known identity 

1 r+oa+ie -, 

— / e" k *e iku dk = g t (u) := -== exp { -u 2 /At\ , e > 0, 

2tt J-oo+ie V4vrf 

we immediately get the following corollary of Proposition |3.3l 

Corollary 3.4. Assume that the maximal isotropic subspace M. satisfies any of the equiv- 
alent conditions of Proposition \2^4\ and defines local boundary conditions on the graph Q. 
Assume, in addition, that the graph Q has no tadpoles. Then the heat kernel of—A(M,a) 
has the absolutely convergent expansion 

\pt(x,y,M,a)] j}j < = Sjjigtixj - y s ) 
(3.12) + W^(w)ft(dist(a i ,«_(w)) + |w|+diHt(y J /,i; + (w))) ) 



The series converges uniformly in x,y £ x 

jesux 
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Proof. It remains to prove that the series in ( 13.121 ) converges uniformly. This follows from 
the estimate 

W M (yv) 9t{dist(xj ,v_(w)) + |w| +dist(x J v,«+(w))) 
<^L= J2 exp{-|w| 2 /4t} 



^7T7 E E exp{-|w| 2 /40 



E rr 1 " 1 ' , exp{-N 2 aLn/4*} 



oo 

< -= £ \I\ n exp{-n 2 a 2 m j4t} < oo, 

v n— 

where a m j n = min{a,}. □ 

In the particular case of a connected graph with X = and standard boundary condi- 
tions, we observe that VUy consists of precisely one walk. Hence, from ( 13.12t we get the 
representation (7.1) in ll28l . which has first been derived in lfl4l by different methods. In 
the particular case Ai = Ai st for compact graphs Q a representation similar to ( 13.121) has 
been obtained by Roth in l47l . 

We will now look at the situation with standard boundary conditions at all vertices in 
more detail. For a given walk w we set 

A?refl(w) = number of times the walk w is reflected, 

-^trans(w) = number of times the walk w is transmitted, 

such that 

Nrefl(w) + iV trans (w) = |w| comb + 1. 

From Corollary [33] we obtain 

Corollary 3.5. Assume that the graph Q is k-regular, that is, deg(u) = k for all v G V, 
and has no tadpoles. Then for standard boundary conditions at each of the vertices the 
heat kernel of —A(M. st ,a) has the absolutely convergent expansion 

\pt{x,y\M s t,a)]jj> =5jjig t {xj -yj) 



(3.13) + E 



2 _ k \ JV r efl(w) /r)X JV tl „ ns (w) 



k J \k / 

■ g t (dist(xj,v-(w)) + |w| + dist(j/ 3 -/, v+(w))). 

4. The Trace Formula 

On the exterior Q cxt = (dV, 0, £, d\s) of the graph Q — (V,T, £, d) we consider the 
Laplace operators A + := A(As = 0, Bs — I) corresponding to Neumann boundary con- 
ditions and A_ := A(As = I, Bs — 0) corresponding to Dirichlet boundary conditions. 

Let J : Hs — > H be the embedding operator defined for any \ £ by Jx — x © 0, 
where the orthogonal sum is taken with respect to the decomposition fi = Tie © Tii, 
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such that JlJ is the identity on fis and J the orthogonal projection in H onto Tig. If 
£ = 0, we set J = 0. 

Theorem 4.1. Assume fnaf f/ie graph Q has no tadpoles. Let the maximal isotropic sub- 
space M. satisfy any of the equivalent conditions of Proposition \2.4\ and assume that it 
defines local boundary conditions on Q. Then 



c tA(M.a) _ Jc tA ±J t\ = L + 1 trjc &{M) T M 

' 2 V nt 4 4 



+ ^ E E^M'klexp {-?»}, t >0, 

where L :— 2^2 jel a i ' s frte fotaZ metric length of the interior of the graph Q, and Wm{i) 
is the weight Wm (w) associated with any walk in the cycle c. In particular, if the maximal 
isotropic subspace M corresponds to a magnetic perturbation of the standard boundary 
conditions (see Example \2.8h then 



e tA(M,a) 



(4.1) 



2v/7rt 




V cG£ pri mPGN V 



where $(c) /s the magnetic flux through the cycle c defined in ( 14.61 l below. 

Remark 4.2. 5i'nce $(c rcv ) = -$(c) ant/ W0w 8t (Crev) = W Mst (c) f^ee (l4TTb below), the 
factor e ip *M m (|4~B can fee replaced by cos(p$(c)). 



Before we turn to the proof of Theorem 14. II we will briefly discuss the trace formula 
(14. U . The first term on its r.h.s. is a familiar Weyl term. In complete analogy with small 
t expansion of the trace of heat semigroups on smooth two-dimensional Riemannian man- 
ifolds [39], the second term depends solely on the topology of the graph: The number 
\T\ — \V\ is the Euler characteristic of the graph viewed as a simplicial complex. In the 
context of metric graphs the Euler characteristics has been discussed in ll27l and [36]. If 
we interpret the quantity | deg(u) — 1 as the local curvature at the vertex v & V, then ( 12. Il l 
gives a discrete version of the GauB-Bonnet theorem for compact graphs: 



E 



ideg^)-!) =\1\-\V\. 



We emphasize that the local curvature at the vertices of the graph is not a curvature in the 
sense of Regge calculus B31 . Regge calculus, however, can be used to define local curva- 
tures on piecewise fiat (or piecewise linear) spaces including Lipschitz-Killing curvatures 
and boundary curvatures [8[. In particular, these curvatures have been used in [8] to give 
an alternative proof of the Chern-GauB-Bonnet theorem for compact closed Riemannian 
manifolds. 

In a similar vein the term \T\ + \£\ — \V\ appearing in (14. U can be interpreted as the 
relative Euler characteristic (cf. |44 |) of the graph Q whenever Q is noncompact, that is, 
when £ ^ 0. In the context of exterior domains in M. d , the relation of Laplace operators on 
forms with absolute and relative boundary conditions (analogs of Dirichlet and Neumann 
boundary conditions) to absolute and relative Euler characteristics, respectively, has been 
established in Q as a relative index theorem in the spirit of ifTHl . 
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The sum over primitive cycles of the graph Q in the r.h.s. of ( 14. U is an analog of the sum 
over primitive periodic geodesies on the manifold in the celebrated Selberg trace formula 
ill (see also JH, EOl, E2, El, EU). 

The remainder of this section is devoted to the proof of Theorem l4.ll 

The heat semigroups e tA± are integral operators with kernels 

hf{xj,yj) ■■= g t (xj -yj)± g t (xj + 

respectively. That the difference e tA ( M, ~> — Je tA± J^ is trace class follows from the fact 
that A(J\4,a) is a finite rank perturbation of A±. For any trace class operator K on Tt 

tr n K= J2 tr-HjPjKPj, 
where Pj is the orthogonal projection in TL onto 7ij. Observe that Pj (e tA ^ M '^ — e* A± ) Pj 

o 

are integral operators on L 2 (Ij ) with kernels jointly continuous in Xj , yj £ Ij (due to the 
uniform convergence of the series in (13.12b ). Therefore, by Corollary III. 10.2 in ifTTl . 
the trace of Pj (e tA (M,a) _ e tA ±) Pj equals the integral of its kernel over the diagonal. 
Hence, from Corollary 13. 41 we get 



Yl / 9t(0)dxj+^2 [9t(G)-ht{xj,Xj)]dxj 
+ ^2 W(w) / 5 t (dist(xj,u_(w)) + |w| + dist(xj,v+(w)))dxj, 

where the sum converges absolutely. We will evaluate the different contributions to the 
r.h.s. of d4.2| > separately. Essentially we will follow the original ideas of Roth developed in 
El- 

1. We start with the terms in ( 14.21 ) not associated with any walk on the graph Q. Simple 
calculations yield 

/ 3 g t (0)dxj = -J= if jel 

Jo V47T< 

and 

/ [fft(O) - hf(xj,Xj)] dxj —T g t (2xj)dxj = =p- if j e £. 
Jo Jo 4 

Summing over all edges j 6 1 U £ we get the following contribution to ( 14.2b 

2. Next we study the contributions from properly closed walks. Let w be a (nontrivial) 
properly closed walk from j € T to j £ X. In this case w_(w) ^ v + (w) and, therefore, 
we have 

dist(xj, u_(w)) + dist(;Ej, w+(w)) = a,j. 

Therefore, 

^2 WOw(w) J g t (\w\ +aj)dxj = ajW M (w)gt{\w\ + a 3) 

properly closed w £ VVj , j ^ 
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Summing over all walks in the cycle c = c(w) and using ( 13.81 ) we get 

Obviously, a j{w) = |c'| if c = pc f for somep £ N and some primitive cycle c'. Thus, 
the sum of the contributions in (I4.2I > from all properly closed walks equals 

Obviously, the relation Wm (pc) = Wm ( c ) p holds for all peN. 

3. We turn to the contributions which are not coming from properly closed walks. In 
this case u_ (w) = v+ (w) and, therefore, we have 



dist(xj, u_(w)) + dist(xj, w+(w)) = 



2^ if «_(w)efl-Cj), 

^•-x,) if «_(w)efl+(j). 
We will call a not properly closed walk w a walk of type A if it is of the form 

w = {jp)Vp,Jp-i,Up-i) ■ ■ ■ , JOjWO) Jo, ■ ■ • ,v p -i,jp-i, v p ,j p }. 

Otherwise a not properly closed walk w is called a walk of type B. By Wf rj and Wfj we 
denote the set of all walks from j to j of type A and B, respectively. 

Obviously, any not properly closed walk is either of type A or B. Any walk of type A 
is invariant with respect to reversion, that is, w rev = w, whereas walks of type B are not. 

The following two lemmas complete the proof of the first part of Theorem l4.ll 

Lemma 4.3. 



jexus W £W^ 



Lemma 4.4. 



iflu£wew fi . Jl i 



Proof of Lemma \4J\ For an arbitrary vertex vo <E V and arbitrary p £ No we set 

d G v (p) { walks °f the form {j P ,v p ,j p -i, ■ ■ ■ , jo,v , j , . . . ,j p -i,v p ,j p } 

with j p e 2}, 

dG £ VQ (p) :={ walks of the form { j p , v p , j p _i , . . . , jo , v , j , . . . , j p -i , v p , j p } 

with j p e £ }, 

5G„ (p) := dGl(p) U 0G* (p), and G % (p) := |J dG VQ (q). 



g=0 



We claim that 



weG„ ( P ) 
(4.3) 1 

- j Y Wm ( w ) erfc 



4 weaG5 o ( P ) ^ 2 ^ 



|w| + 2oj (w ) 



HEAT KERNELS ON METRIC GRAPHS 



17 



holds for all p e No. Here for brevity we set & Vo := &(A4 Vo ), where M Vo C d C Vo is the 
maximal isotropic subspace from the orthogonal decomposition (12. 1 9b . 



erfc(s) 



e _u du 



denotes the complementary error function JT|. The proof is by induction. For p = we 
have 



Wm(w) / 9t{2xj(w) + |w|)ttej( w ) = - tr Ct , o [& V0 ] 
\ £ ^(w)erfc(^). 



WSG„ (0) 



wG9GJ Q (0) 

Now assume that ( 14.31 ) holds for some peN and consider 



( w ) / 9t(^ x j(vr) + IwD^j^w) 
vGG„ (p+l) Jl n*>) 

Y Wa^(w) / St(2a; 3 -( W ) + |w|)^ i(w ) 

S„„(p+1) ,// J(w) 



wSG„ (p) 



veaG„ ( P +i) 



it«,[ej-i E ^(w)e rf c(W±|i«) 

i £ ^(w)erfc(W) 



rfc 



- Y Wm(w) 

we9G5 (p+i) 

|tr £ „ o [6, ]-J ]T ^(w)erfc 
weaGj ( P ) 

i £ ^(w)erfcfW) 

wG9G„ (p+1) 

- ^ W M (w)erfc 



2a 



j(w) 



2 V* 



Iwl + 2a 



'i(w) 



2Vt 



wl + 2a 



weOGJ o ( P +i) 



i(w) 



2Vt 



It remains to prove that the sum of the second and third terms on the r.h.s. is zero. Let w S 
dG Vo (p + 1 ) be arbitrary. Write the walk w as w = { j p+ 1 , v p+ 1 , j p , . . . , j , v , j , . . . , j p , 
v p +i,jp+i}- Then w' := {j p , jo,v ,j , ...,j p }€ dG Vo (p) with j p £ I. Hence, 

WmW = [&v p+l } jp , jp+1 [e Vp+1 } jp+1 , jp w M (^') 
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and |w| = |w'| + 2aj( W '). Thus 



J2 T4^(w)crfc(M) 

/eOG„ ( P +i) vzvl/ 



w'| + 2a,-( w /) 



w'e9GJ (p) 

jp+ie5(u p+ i) 

By (iii) in Proposition ^. 41 we have &1 p+l = I and, therefore, 

[®fp+i]i P o'p+i [®fp+i]i P +i jp = [®Dp+Jipjp = 

ip+iS5(t)p + i) 

which completes the proof of (14.3b . 

By the absolute convergence of the series ( 14.21 ). from (14.31 > it follows that 

^So X! w m(w) 5t( 2 ^j(w) + |w|)rfx j(w ) = - tr Cvf! &v - 

P_>D °weG„ (p) "^"(»> 

Observing that 

P ^°° v EV v,<=G V0 (p) Jl n™> 
we obtain the claim of the lemma. □ 



Proof of Lemma \4~4\ Any not properly closed walk of type B is obviously of the form 

{jpiVpJp-l, . . . ,jo,V ,5,V ,j , . . . , jp„l,D p , j p } 

for some p £ No, where 5 stands for the sequence of internal edges and vertices v[, . . . , 
v' n , i' n with %' x ^ i' n . For an arbitrary p £ No we set 

dF s,v (p) '■ = {walks of the form {j p , v p , j p - X , . . . , j , v , S, v , jo , . . . , jp-i . % , j P ] 

with jp £ I}, 

dF s,v (p) : = {walks of the form{j p , v p , j p -%, . . . , jo,v ,S, v ,jo, ■ ■ ■ , j P -i, v p , j p } 

with j p ££}, 

p 

dF s , V0 (p) := 8Fl vo (p) U 3F*„ (p), and F, fCo (p) := (J <?F S ,„ (g). 

9=0 

We claim that 

V^m(w) / g t (2x i(w) + |w|)dx i(w ) 
1 v-^ , s / Iwl + 2a, 

3 „§ , 1 M '- M(w) ( 27T 



(4.4) 

x j'( w ) 
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The proof is again by induction. For p = we have 



( w ) / St(2£j(w) + IwDcfe^v 
■ w. ...mi u ' / j(w) 



weF 3 ,„ (o) 

(4.5) = w m{™) / fft(2x j(w) + |w|)dx i(w) 

+ X! w m{w) I gtfajfr) + |w|)cfcEj-( w ). 
weap£„ (o) "'^(w) 

Obviously, for any w g i*i lt , o (0) we have WOw(w) = [6^]^ [St) ]i^,jWjU' where W' M 
is a weight associated with the sequence s and j — j(w). Therefore, if j G £, then 



and, if j £ I, then 



= llG^i&v^jW^ (erfc (^\ - erfc 



|w| + 2 % 



2Vt 

Again we use &1 = I, which in combination with i[ ^ i' n gives 



3&S(v ) 
j65(«o) 

Combining this with (14. 5t , we get the claim (14.4t for p = 0. The proof of the induction 
step follows the same line as in the proof of Lemma |4~3l and will, therefore, be omitted. 
By the absolute convergence of the series ( 14. 2K from J4.41 i it follows that 

^ X! W m(™) / 9t(2x ]{vr) + \vf\)dx ]{w) =0, 
P ^°°weiWri 3(w) 
which completes the proof of the lemma. □ 



To complete the proof of Theorem 14.11 it remains to consider the particular case of 
magnetic perturbations of standard boundary conditions (see Examples 12.61 and I2.81 l. We 
assume that the maximal isotropic subspace M. corresponds to the magnetic perturbation 
of the Laplace operator — A(A4 s t , a) with standard boundary conditions, that is, 

M=M u st with U=®U V , U v = diag ({e^W} ieS(tl) ) . 
vev 

First we calculate tr*; &(M). Using (l2~TT i, ( I2.241 i, and ( 12.251 ) we get 

tr K 6(M) = tr K S(M st ) = ]T H [©«(<Hrt)lw 

= ^(2-deg( V ))=2|y|-|f|-2|2:|. 
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Let Hi (Q, Z) be the first homology group of the interior Q- m t of the graph Q. There is a 
canonical map T : <t — ► H\(Q, Z), which satisfies 

r( Cre „) = -r(c), r(pc)= P r(c), P eN. 

In particular, T(c) = when c = c rcv . Therefore, the map T is not injective. In general it 
is also not surjective. For any cycle c G £ we set 

*(c) := *(T(c)), 

where $(F(c)) is the magnetic flux through the homological cycle T(c) as defined in ll27ll . 
If 

W = {j ,V 0> ji, ■ ■ -Jn^nJo] 
is an arbitrary walk in the equivalence class c, then by explicit calculations the magnetic 
flux through the cycle c (see l27l ) is given by 

n-l 

(4.6) $(c) = (VjuiVk) ~ ¥>j h+1 (Vkj) + {<Pj„(v n ) - <Pjo( v n)) • 
Obviously, 

ff M (c) = ^(c)e i$ W 

and 

(4.7) $(c rev ) = -*(c). 



By Proposition |231 we have Wm(Zi-cv) — Wm{i)- In particular, if c rcv = c, then W^(c) 
is real. Since $(c rev ) = -$(c) and W Mst (c rcv ) = W Mat (c), the factor e ip0 ^ in ( H3T > 
can be replaced by cos(p$(c)). This completes the proof of Theorem l4.ll 

5. Applications to Inverse Problems 

In this section we present an application of the trace formula in Theorem |4.1| to inverse 
spectral and scattering problems. Throughout the whole section we will assume that the 
maximal isotropic subspace M. satisfies any of the equivalent assumption of Proposition 



For the noncompact graph Q let S(X;A4,a) : JCg — > K,£, A > 0, be the scattering 
matrix for the triple (— A(A4;a), — A^,^) defined in ll26l according to the scattering 
theory in two Hilbert spaces lISTl Chapter 2]. Here J is the identification operator defined 
in Section |4] The scattering matrix is continuous with respect to the spectral parameter 
A > (see J24) or J29] Theorem 3.12]). 

Let £(A;.M,a) be the spectral shift function associated with the triple (— A(vW;a), 
— A + , J) (see [51 , Section 8.11]). It satisfies the trace formula 



tr w 

(5.1) 



e tA(M;a) _ J c t&+ jt 



tr H£ [(jij-I H£ )e tA +] 

- -t / £(A;.M,a)e~ tA dA, t > 0, 
Jo 



and is fixed uniquely by the condition £(— 1; M, a) = 0. From the definition of the opera- 
tor J it follows that the second term on the r.h.s. of (I5.lt vanishes. Thus, 

tr„ (e tA ( M '*> - Je tA +jA = -t [°° e" At ^(A; M)d\, t > 0. 



By the Birman-Krein theorem the spectral shift function is related to the scattering matrix, 
(5.2) det K : £ 5 , (A;A^,a) = cxp{-2ivi£(\; M , a)} a.e. A e K+. 
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By the continuity of the scattering matrix and due to (15. 2i one can choose the branch of 
the logarithm such that 

(5.3) s(X;M,a) := ^ logdet K ; £ S , (A; M, a) 

is continuous with respect to A G (0, oo) and satisfies 

s{0+:M,a) = -it (£,(0+:M,a) + N(0+;M,a)) , 

where N(X; JA, a) is the counting function for the eigenvalues of the operator — A(A4 , a). 
The function ( 15.31 ) is called the scattering phase. The scattering phase and the eigenvalue 
counting function uniquely determine the spectral shift function, 

£(X;M,a) = --s(X:M,a) - N(X;M,a), X G R+. 

For compact graphs the spectral shift function is determined by the eigenvalue counting 
function alone, £(A; M,a) = -N(X;M,a), X G K + . 

Proposition 5.1. Assume that the graph Q {compact or noncompact) has no tadpoles. Let 
the maximal isotropic subspace M satisfy any of the equivalent conditions of Proposition 
\2.4\ and define local boundary conditions. Then the spectral shift function M+ 3 X i— > £(A) 
uniquely determines the set 

(5.4) [t >o| Yl W M (c)¥>0}. 

cGC 
\c\=t 

In 1 35 1 the set ( 15.4b is called the "reduced length spectrum". 

Proof. Using standard formulas for the inverse Laplace transform and, in particular, the 
fact that t~ z / 2 e~ a / 1 , a > 0, is the Laplace transform of {■ko)^ 1 / 2 sin(2VoA), from Theo- 
rem l4.1l we get 

-£{*M,a) = -VX+ ]tr K 6(M) - ^ 

7T 4 4 

+ ^ E 51 -W^(c)Psin (p|c|VA) , A>0, 

where the series converges in the sense of distributions on M + . For k G K define the 
function u(k) via 

'-£(k 2 ;.M,a) if k > 0, 
Z(k 2 ;M,a) if k < 0. 
Using ( 15.51 ) we can calculate the distributional derivative of u, 



u(k) := 



«'(k) = ^+Qtr K 6(M)-H) S(k) 

+ ~ E E^W^IclcosHclk), kG 



7T 

CSCprim PGN 



where 5 stands for the Dirac ^-distribution. Its Fourier transform with respect to k yields 

J e iujk u'{k)dk = 2L5(u) + i tx K &(M) - ^ 

+ E Wa4(c)"]c| - p|c|) + tf(a; +p|c|)] , 
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which implies the claim. 



□ 



Theorem 5.2. Assume that the graph Q ( compact or noncompact) has no tadpoles. Let the 
maximal isotropic subspace M. satisfy any of the equivalent conditions of Proposition \2.4\ 
and define local boundary conditions on the graph Q. Assume, in addition, that 

(i) the lengths di (i € T) of the internal edges of the graph Q are rationally indepen- 
dent, that is, the equation 



2 , n i a i — 

with integer n, € Z has no non-trivial solution; 
(ii) for any vertex v G V, none of the matrix elements [& v ]j,j'> € S(v), vanishes. 
Then the spectral shift function M + 3 A i— ► £(A) uniquely determines the interior Q- lrA of 
the graph Q. 

Proof. All arguments of Section 4 in [35 1 remain valid for boundary conditions satisfying 
assumption (ii) of the theorem. Thus, the set ( 15 ,4b uniquely determines the interior Q lnt of 
the graph Q. Combining this with Proposition l5.1l we obtain the claim. □ 

Note that assumption (ii) of Theorem l5.2l implies that if A4 V corresponds to the standard 
boundary conditions at the vertex v or its magnetic perturbation, then necessarily deg(i>) 
2. 

Remark 5.3. As in B42II the assumption on the rational independence of the edge lengths 
can be slightly relaxed. 
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